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This part concerns the height-gain differential equation in order to obtain a series for 
the complete solution which starts with the extended W. K. B. approximation discussed in 
part I. The coefficients of this equation depend for each mode amongst other things on the 
parameters A,- fixing the refractive-index profile. However, the explicit dependence on 
these parameters can only be given in terms of expansions with respect to (te)"^^^ (A:oa = cir- ' 
cumference of the earth divided by the wavelength). In turn these expansions are derived 
with the aid of other ones for the complex turning point connected with the height-gain 
differential equation. T\u} final expansion for the solution of the differential equation is 
substituted in the boundary condition at the earth's surface. This leads to corresponding 
expansions, with respect to {koa)~^^^, of the quantity ui{a), and next of the eigenvalues / 
themselves. 

16. Canonical Form of the Height-Gain Differential Equation 

According to section 'A the field of a dipoie situated in a spherically symmetric stratified 
medium can be expressed in terms of a single scalar quantity 11. The latter can be expanded 
in modes Ji{r)Pi (cos 9) the radial parts of which have to satisfy the height-gain differential 
eq(8), viz, 

{^2+'clrnj{r)y{rJ,(r)}=0. 
According to (9), (14), and (15) rui is defined as follows hi terms of the profile ^eff=(^/^)^eff' 



m,(r)=^VM^ff(r)-C'f. (65) 

Moreover, the height-gain functions should fulfil the radiation condition at nifinity, as well as 
the boundary condition (23) at the earth's surface. 

Following the methods of Langer [1937] and Pekeris [1946] we hitroduce the new variables: 



and 



Ui{r)=ko r mi{s)ds=ha {' -^Ml^{s)-C] (66) 

Jri Jri ^ 

G>=^/^^)Mr)=J^ [Mlnir)-C\Y'\ft{r), (67) 

tt If ^ 



the former of which has been discussed in section 5. The influence of the refractive-index 
profile enters explicitly by the occurrence of Meff, and implicitly by that of the turning-point 
level Vi defined by (14). The lu)ight-gain differential equation then proves to be equivalent to 
the new equation : 

^+P,(«0(?,=0, (68) 
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in which 

4 mf 2 Ml 



^^(^0-1+7^-F-^-^- (69) 



The eq(68) will be the starting point for our further investigations. 

17. First Expansion for m.i{r) as a Function of Ui{r), With Coefficients Depending 

on Ti 

In order to solve (68) we have to determine Pi(ui) as a function of Ui. This can be done 
for each eigenvalue I in terms of an expansion the^coefficients of which depend on the parameters 
A^ of (17). The latter characterize the refractive-index proj&le. In view of (69) we start with 
the determination of m^ as a function of Ui. The tedious derivation of the coefficients of the 
corresponding expansion will be given in steps, to be discussed in the next sections. 

We first notice that m? could be expanded into a Taylor series in r— r^, since mj is a contin- 
uous function of r according to (65). However, this series starts with the linear term; in fact, 
the relation M^ttiTi) = Gi (see (14)) involves mi{ri) = 0. Weinfer a proportionality of mi(r) with 
(r—ri)2 for small values of r—vi. In view of (66) this involves a corresponding proportion- 
ality oiuiwith {r—riY^^j or of r—r^ with t^^^^, and therefore of m? with ^^5/^ Similar considerations 
show the proportionality of the next terms of the expansion of mi (in terms of Ui) with t^?/^, 
u^i^^,uy^, etc. We therefore put: 

5=^om,=lg5-L_, (70) 

the dimensionless coefficients B^l^ of which are to be determined next. In this section these 
coefficients will be connected with the properties of the profile near the turning point Vi. 
We first multiply the expansion (70) by itself, so as to obtain 

We substitute s = m—j, and invert the order of summation with respect to j and m. The 
resulting expansion reads : 

CO ^yl(Wi+l) m 

kWrnl=T>7T-^-=Tyi:BpB»lr (71) 

m=o wci;^'- ^ j=0 

The coefficients Bj^^ can now be obtained from the derivatives of (71) with respect tor, 
taking r=ri. For instance, the first order derivative yields, applying both (70) and (71), 






Since r=ri involves Ui=0, we have only to derive the term m=n=0 in order to get the value 
of the left-hand side at r=ri. Similarly, a further differentiation of (72) with respect to u, 
multiplied by (70), yields the second derivative of mj; its value at r=ri reduces to the constant 
term of this new expansion, and so on. 

We mention the expressions thus obtained for the first few derivatives of m^ at r=ri: 
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From these relations we may solve in succession the coefficients 5*1,' , B^[\ etc. We obtain : 

^^ -56\3r7;/ Ir?; 2577;^/' 

5^^>-/^"V"^|-^-^^^^+-^^\ (74) 

^"1377;/ \l6r7; 7.?;^ +2800 r;;^ / ^^^^ 

18. Second Expansion for mi{r) as a Function of u/r) With Coefficents Depending 

on the Profile-Function M{r) 

According to (73) and (74) tlie coefficients B^^ of (70) depend on the position of the turning 
point Vi which mnj l)e derived from (18). In this latter expansion the parameters Ay are param- 
eters known from the given profile of the refractive index. In order to obtain the turning 
point from (18) in its complete explicit dependence from the profile, we shall have to Imow the 
coeflficients mi{a) entering in (18). Therefore, we next want to obtain a connection between 
the latter and the profile i)arameters A^. It can be arrived at from an expansion of m.i{r) in 
terms of i^z(r), to be derived as follows. 

We start from the expansion (21). In view of (65) it can be represented by 

Ui{r)=ha iz af;(r)m,(r)2^■+^ (75) 

;=0 

in which 



a.W-(-l)^jff ...2j.^3^^^f(,)j,^ ^cZMIg y. (7(>) 



[M- 



We invert this series, taldng into account that the first term proportional to m\ in the original 
series involves a first term proportional to u^^ in the new one. The complete inversion series 
can be put in the form : 



»..«=«.{^}'%^',<-){=|^}"%c,,,{^}"V..«{^} 



7/3 

+ ... (77) 



In fact, substituting this expansion itself ajid the corresponding expansions for its suc- 
cessive powers, such as 

(78) 



-'W=^o^{£}^%2C„r>,{|-;}^%(6?+2r„Q{^;}^%2(c„6'3+6^,Q{ii 



into (75), we find by equating the corresponding powers of Ui in both members unique relations 
for the determination of the new coefficients Cj{r) from the other coefficients aj{r). The first 
relations of this type read evolicitly: 
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c\{ry. 



[aoir)}' 






C,{r)- 



4al(r) — 3aoir)a2ir) 



-65 



9a'o"'(r) 



Csir)- 



alir) 



-90: 



«^(r) 



-27 



«2W 



(79) 



These coefficients thus determine the expansion (77) for /«;(/') in terms of Ui{r). 



19. An Expression for the Position of the Turning Point in Terms of the Param- 
eters u,{a) and Ay 

Sucli an expansion can now be derived by applying (77) to r=a. We find 

^.^ fu,(a)Y'\ ru,ia)Y'\ fu,(a)y'\ /«,(«) \^'^ ,„„. 

^'(«^=^°ii^| +^'\^^/ +''^1^^/ +^^w/ •■•' ^''^ 

where yj = Cj(a) follows from a substitution of r—a in (79). Introducing the new quantities 
fij=aj{a) we then have: 



7i = 



72= 



73 = 



3^r 

_R(^ ^4- on ^—17 ^ 



(81) 



81^J« 
Further, these coefficients 0j are connected as follows with the profile, in view of (76) for r=a, 



'^' ^ ■'357 (2j + 3) 



rf^ 



1 



c?{Myr-)}'-^^ilf|g 



(82) 



All derivatives occurring here can be expressed in terms of higher-order derivatives of AP. 
For instance, we find : 

d 



1 

dMi^ dML 



1 



1 dr' 



dr 



\,/dm«y r/dMln^' 



dmi. 



fdmiA ' 

_ 3 \ rfr^ / 

{dMl^ymi, ^,/dM!,V'r^/dM!,V r/dMinV' r /clM^' 
dr \ dr J \ dr J \ dr / \ dr ) 



3 dr"" 1 dr" 
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We only lu'cd all tlicsc quantities at r=a in order to evaluate the coefRcicnts /S^ from (82). 
All terms then occunijig; can be expressed in terms of the profile constants A^-, by applying the 
defijiition (17). Tlic first few prove to be given by: 

2 



8 / 2 3A2-A3 3A1\ 



. _ 16 / 6 , 12A,-4A3+A4 , ISAj-lOA.Aa , 15A^,\ 

'^='-945 Vrr A| + i\ +^/ ^^'^^ 

The coefficients 7^ of the expansion (80) can tlius be evaluated juunericallv by first computing 
tlie ^/^ with the aid of (83) from the profile constants A^, and next applying (81). For a further 
evaluation of the expansion (18) for tlie position Vi of tlu> tur.ning point, we further Juive to 
determijie the even powers of the expajision (80) in order to arj-ive at iui exphc^it i-epr(>sentation 
of ri in terms of both A^ and Ui(a). Obviously we then get aji expansion of {ri—a)/a hi terms 
of increasuig integral powers of th(^ cpiantity 



{ Ui{a) \ 
ha J 



2/3 



with coofficieJits depending explicitly on both jj and Aj. Tliis rather tedious procedure leads 
to the following final representation: 



n— 



a fuiia)^' fui{a)Y^' Cu,{a)Y^-\ fuiia)^^' ,_., 

~='\-ha-f +^ni^/ ^''Vf^f +''Vha~f •••' ^^^) 



a 
the new coefficiojits of which can be derived as folloAvs from the coiistajits determhied previously: 

^'-'-~~^~2Af^»' 

7i+27o72 2A2 3 /A3 Al\ , 

P'= A^ A[^«^'+(6Ar2Xf>»' 

c 7073+7172 A2 ,„ 2 210 3 \ 1 /^3 3Ai\ 5 , / 5 A2A3 5 Al Aj \ s , , 
P.= -2 —j-^ ^ (37o7f+27o72)+(^-^J7o7i+(^Y2-Ar~8 Ar24A!J "">■ ^^^^ 

Summarizing, all coefficients of the expansion (84) can be evaluated numerically by first 
deriving the ^/s from the profile parameters Aj with the aid of (83), next tliCT^'s from the jS/s 
with the aid of (81), and then the p/s from the 7;'s with the aid of (S5). 

20. Successive Developments Resulting in an Expansion of the Coefficients 
Bn{l) in the Expansion (70) for the Dependence of m/(r) on Ui{r) 

Such an expansion is required for the subsequent analysis enabling tlie numerical evalua- 
tion of the eigenvalues and the corresponding expression of the field components. We start 
with a derivation of a series for the parameters r?/*^ of (73) hi terms of tlie turningpoint co- 
ordhiate Vi, 

Let us first differentiate the square of (65). The definition of r;J then uivolves 
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which reduces, applying (14), to : 
The Taylor expansion : 



MUT)=±r^,{r-ay (86) 



may be used here for a further reduction which results in : 



rf^oa^Xi-l)! 



The higher-order coefficients »?;*' are obtained in a similar way. For instance, by differ- 
entiating twice the square of (65) we obtain first: 

.;'-^£g(Myr)-c^}]^^^^ 

which can be reduced with the aid of (14) to: 

When substituting the first and second derivative of (86) we arrive at: 

„_ V 4 ^ A, (n-g)^--' I 1 f. A, {r^-ay-n . 
^' ""L T\n<^' ij-iy. ^rjika' (i-2)!j" 

By a continuation of this procedure we get expansions for all coefficients 77^^^^ in terms of 
powers of ri—a, divided by another power of r^. In its turn, the latter power can be developed 
with the aid of the binominal expansion for 

rf a^ \ a J 

When assembling identical powers of ri—a we thus obtain expansions of the form: 

Vi'''-±qj'^y> (87) 

the lower-order coefficients of which prove to be given by: 



A, 
2io=Ai; gii=A2— 2Ai; 2i2=3Ai— 2A2-f-Y; 



gi3=-4Ai+3A2-A3+|^; 



2i4=5Ai-4A2-h2 A3— 3-'+2|; 
220=— 4A1+A2; g2i=12Ai— 6A2-I-A3; 222= — 24Ai + 15A2— 4A3-l-y; 

223=40Ai-28A2+9A3-| A4+^; 
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224=-60A,+45A2-16A3+^A,-|-=+^«; 
230= I8A1-6A2+A3; g'3i=-72Ai+36A2-8A3+A4; 

232=180Ai-108A2+30A3-5A4+^; 

g-33=-360Ai+240A2-76A3+15A4-2A5+^'; 

234=630Ai-450A2+155A3-34A4+^ A,-j^ A,+^' 

24o=-96A,+36A2-8A3+A4; 

g:4i=480Ai-240A2+60A3-10A4+A5; 

242=-1440A,+840A2-240A3+45A4-6A5+y; 

7 A 

243=3360Ai-2160A2+680A3-140A4+21A5-~A6+-^'5 

o u 

g44=-6720Ai+4620A2-1560A3+345A4-56A5+7A6— |a7+^- (88) 

As a next step we may substitute in (87) the expansion (84) and the successive powers of 
this expansion. A new power series for jji*^, viz 






then results. Its lower-order coefficients can be represented as follows: 

tnZ=g.nlP3 + 2qn2PlP2 + qn3P\', 

tni=qnlPi + qn2{pl + 2piPs)+^qn3PlP2 + gn^pl (89) 

If these coefficients are also known we are enabled to derive also expansions of the co- 
efficients Bi^\ defined by (74), in terms of the quantity Ui{a)Jkoa, For instance, we have 

-(If VTi , ltnSn,{a)Y'^ fit,, ltL\fu,{a)\''' 1 



and 






4/3 



10 ij 10 

t,o+t 



_3 <2or I r^2i tn\ futiaU'" 

lO^.oL \<2o iioj \ ^oa / 
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and so on. In this way we get in general expansions of the form; 

The first few^ of these new coefficients are connected as follows with the former quantities tnj'. 

7^ 3 J ^22 '^11^21 ^20^12 I ^20^11 \ • 

iU ^fio fjo flO ^10 J 

T) _3_ J t^Z ^11^22 ^12^21 '^13^20 I r> ^11^12^20 | ^21^11 ^20i^ll \ . 

^13— io"S / " /2 "~ /2 "~ /2 "T^ 73 T /3 "" /4 f ' 

'" 56V3<,o/ \<,o 25«fo/' 

7? ^ 9 / 2 y r<3i 4 tsptu 38 <2o^2i , 133 tyn\ . 
'' 56V3^,o/ \<io 3 tio 25 <?o"^75 f?„ /' 



*^12^30 4 ^3i/i] , t 



+f 

f 10 



I 38 ^11^20^21 38 ^20^11 19 V tio/ 

75 610 75 i^io 25 /io 

I 38 ^20 / ^i2^2o ^11^20 ■ tni2i . \ 

^5 610 \ C 10 ^10 ^10 / 

\9 ^10 3^10/ V^io 25 tiQ/ J 

JD _./ 2 \ J i^40 ^20^30 I 213 ^20 1^ ' 

'" \3<io/ \l6<,o 7<?o 2800<io/' 



5 ^40^11 ^20^31 ^21^30 I §_ fiit2ots0 

10 48 ^10 7^10 '^^lo 21 ^10 



639 1^20^21 781 ^11^20 I /Q-|\ 

"2800 tlo 2800 -4 ^* ^ -^ 



^11^20 \ 
t\o J 



These coefficients, occurring in the expansion (90), determine the dependence of mz on Ui accord- 
ing to (70). 
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2 1 . Explicit Form of the Differential Equation (68) 

Tlu^ expansion (70) for m^ as a function of Ui can also be differentiated in order to obtain 
corresponding expansions for druildui and d^mildu\. It is also possible toj deduce from the 
latter developments of the quantities 

^2 d^Mi 



/ dmi \ 
\dui) 



1 dui 
- and 



mf m 



occurring in (69). This elementary formal procedure leads to the following expansion for the 
coefficient Pi of the final differential equation (68): 



^^(^) = 1+T^+Z; • — ^ — ui^^-'\ (92) 



Its first few coefficients read: 



Co"' 



BP' 



^2 



"9 B",''' 



C[''=-, 



SBi'^B^'' ABi'' 2B[ 



ar 



9 B^'^' 3 5^^^ 9^(03 
___8 5i^ , 8 Bi'^" , 2B{'^Bi'' 14 Bi'^'B^' 



^,j^_40 5|^ 28 B^B^ 32 5{[^_8 ^B{^\^_26 B^^B^ 22 Bj^'^Bi'' 4 5{^ 
^^ ~ 9 i?i^^+9 5^^ +9 i^m^ 3 ^(0^ 9 5^' 9~5r' 9 j?^^)^^' ^^"^^ 

when derived straightforwardly with the aid of (70). In its turn the quantities j5^" are 
connected with the eigenvalue parameter Ui(a) according to (90). The later analysis requires 
the explicit dependence of the coefficients Ci^^ on this latter parameter. This dependence is 
obtained by determining the powers of the expansion (90) and also the relevant fractions 
resulting from the divisions of such powers. 

The evaluation of this procedure leads to further expansions of the form: 

the new coefficients of which may be computed by performing the mentioned procedure nu- 
merically, without applying the complicated analytical expressions which connect the Ck/s with 
the coefficients Bnj determined previously. However, for convenience we give the expressions 
for the simplest coefficients Cjcji 

n — 4 ^20 , 1 ^10 



9 5„, 



-"oo 



^_ 4:Bn,4:B,oBo^+2BioBn 2 B,,Bi, 



^''~~9B 



9£„„^ 9Blo 9 Bio 

4B2, , i{B.nB,, + B.,oBo2) + 2BioB,2+Bj, 4:B,oBl,+B,o(4:B,nBn+2B,oBo,) , -Bfo^g, 



4530, 85io52o_25|o 

«>- 3B0J9 Bl, 9 Bio' 
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3 Boo Bio Bio 3 B, 



^r~f 



2BioBso+-B2o 

^20— o I? "I D2 n Z?3 ^o 7?4 V^^/ 



8^40, ^"^"^9 ^^ 14 ^^0^20 . 1 5to 
■3 5oo"^ Bl, 9 5^0 



22. Solution of the Differential Equation (68) 

This equation depends completely on the function Pi which is represented by (92) as a 
power series in {koa)~'^'^. The only term remaining for koa-^^ involves the solution G{u) 
='yJuHi]l(u) known from the ^'extended W.K.B. approximation.'' The next terms may be 
considered as corrections for large finite koa and suggest the following representation for the 
complete solution: 

We substitute this series into (68) and assemble the contributions proportional to one and the 
same power of (AroCt)"^''. It leads to the following set of equations: 



du' 

d% 






m-2 



lS+('+i?) ^»=S ^'^-'^-^-''"''—''OM. («>2). 

Starting from the normalized ^^extended W.K.B. approximation'' go{u) =^fuH[%{u) , satis- 
fying the first equation, we may take ^i=0. Taking into account these special values of ^o 
and ^1, we find the next equations: 

etc. All these equations can be solved easily in succession, by substituting for the unknown 
function in any further equation the sum of (a) the Hankel function H[)z{u) multiplied by a 
finite number of (positive or negative) powers of v}'^ with indefinite constant factors, (b) the 
Hankel function H^]l(u) multiplied by another number of such powers. The powers to be 
introduced here, as well as the coefficients in front of them, can be determined conveniently 
by applying the following identity: 

40 



The resulting functions g^, substituted into (96), yield the following first terms for the 
complete solution of (68) that starts with the ^'extended W.K.B. approximation": 



(-| Cl+^^ umi%{u)+^ e-'^''C,u"mi%{u) 



8/3 



23. Equation for the Eigenvalues in Terms of the Functions G and nrtj 

This equation results from the boundary condition (23) holding at tlie earth's surface 
(r=a). It reads as follows for the individual modes lLli=Ji{r)Pi{cos 6): 

^(^•/i)=-^/^at r=a. 

In view of the definition (67) for the function G this proves to be equivalent to: 

drrii 

IbG 1 dr r ^ ^^^^. 

y^-^ o =~ at r=a, (98) 

We shall pass from /' to the variable u which can be performed with th(^ relation 

dui J 

resulting from (66). The equation (98) then transforms into: 

IdG. hadmi _ ^ ,.^. 

24. An Explicit Form for the Equation for the Eigenvalues 

The function G{n) is represented explicitly by (97). Its derivative dG/du can also be 
expressed in terms of tlie functions Hi/^ and H2/3, using the identities: 

We tlius find the following expansion corresponding to (97): 

'^2 ^" (haF" '^^ ~ (hay 

{(~i ^'°+^) '*'"+^ C,u'"'yH^^liu)+e-'"' (I <7^+^) u^'^H^^iu) 
+ ^^^^p + . . . 
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A division of (97) by this new expansion leads to the following corresponding representation 
for {d6ldu)iG: 



d6 ^,l/3/l_g.W3^l!W\ 






(A-offl)^ 






8/3 



In order to evaluate the first term of (99) the latter expansion has to be multiplied by that 
of koanii which, in its turn, is obtained (in the special case r—a) with the aid of (70) and (90)- 
We find, for each eigenvalue I 



2/3 



+ (Bo3+-Bi2+-B2i+Z?3o) Ye-^,+ {B,,+B^,+B,, + Bn+B^) ^ . . . for u=u,{a). (100) 

The product constituting the first term of (99) can thus also be expanded with respect to 
(kiiO)'^'^. The corresponding expansion of drui/du, entering in the next term, follows from the 
derivative of (70); it depends on the same coefficients Bnj of (90) since we have to take r=a 
throughout. By these procedures we arrive at the following final form for the equation (99) : 



r = -.--i?o„.-^ (M-+{aoo+«o.e--.||g} 



^H^IM,,.,„_ 



«3o+a3.«^+(e-'-^%32+e'^%33)«ii^l + e-'-^%34ti'Si7S 



(koa) 



{k(,aY 



(101) 



in which «is short for Ui{a) whereas the new numerical coefficients introduced here are defined by: 
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^00 — 7^ ) <^01 — (^Ol + ^io), 






(ln= — {Bo2+Bii-{-B2o), 



^12 — o ^OO^OOj 

«2o=§-" (3Co,+(7.„)+| (i?o,+5,o)Co„-^^-| 5.,-i^.„, 

^21 = -OqoCio, 



023= - / ^00 (I a. +%) + ! ( ^^1 + -K.o) ^00 } ' 

^30= ^00 ( ~~4 Q)0+7T ^20 )^ 

«'31=^Oo( YO ^'^^^"'"2 ^°^~' 2~/~' ^^ ^(3f^01 + C'lo)+:^ (^02 + -Bll + ^2o)C'oo 

-l(B,,+2B22+SB,,+4B,o), 

a32=i>oo ( 7 ^20 T^oo)^ 

BooCn + (Boi + Bio) C'lo 

^33= , 

(i'ZA= — {Boi+Bis+B22+Bsi+B^o), 

as,= -[^00 (^ C^2o+§-'+| a2)+(^oi + 5,o) (I Co,+^^)+| {Bo2 + Bn+B2o) ^oo} 

a36=^5ooC^oV (102) 

The equation (101) fixes the eigenvalues I of our ])r()l)lem since it determines the value of 

u,(a)=koa r^^jMUs)-C^, (103) 

J ri S 



Avhicli leads to a well defined value of (7z=V^(^+l)/^o^- This will be discussed further in sec- 
tion 26. 

25. Solution of the Equation for the Eigenvalues 

A first approximative solution of this equation is obtained when neglecting all terms but 
the first in the right-hand side of (101). The corresponding value % of Ui[a) satisfies the 
equation 
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if we put r = (Jk^af^V . In the case of a homogeneous atmosphere the relation (104) constitutes 
the base for the conventional diffraction theory for the propagation of radio waves round a 
smooth spherical earth. 

The decrease proportional to powers of (^o^)"^^^ of the next terms in the right-hand side of 
(101) suggests to expand the rigorous solution as follows, for each eigenvalue u=Ui{a), 

''"''^+(M^"^(M^''"^(M^ • • • ^ ^ 

It is then necessary to expand all coefficients in (101) accordingly, with the aid of Taylor 
series around u=u^. This will be indicated in detail for the first coefficient, viz, 






After substitution of (105) we get an expansion with respect to (koa)'^^^ that starts with: 






■ 1 (u\ d^ . ^, d\f u'm^%(u) \ 



1 fu\ d' d' , d\f 'u"mi%(u) \ , ,,_, 



All higher-order derivatives occurring here can be expressed in terms of the Hankel func- 
tions of the orders 1/3 and 2/3, by applying a number of times the elementary rules known from 
the theory of Bessel functions. For the evaluation of the present coefficients we need the 
relations 

du H\]\{u) -^ y ^' H[f{u)j' 

du' Hi%{u) -3u"^^^'' Hi%{u) 2u'''Hi]l\u) HDUuf 



<i_ n^'mUu) ( 1 Y 2 m%{u) /-,^V3 , J_\^ 



"in) 



2e*^'^m]t{u) ,^.y,mi\u) 



&u^'^- 



u'l' mtiu) m%\u) . 

In all these coefficients we may substitute, in view of (104), 



H{%{u,) Boo u\' 



In this way we get the following expansion replacing (106): 
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Corresponding expansions are found for the other quantities entering in (101). They start 
with : 
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Substituting all these expansions into (101), and equating to zero the total coefficient of 
one and the same power of {koa)~^^^ [with an exception of the power {k^aY^^ which leads to 
r'=r'] we obtain relations from which Ui, U2j . . . can be derived in succession. In fact, the 
u with highest suffix always occurs linearly so that we only have to solve linear equations in 
order to derive Ui from Uq, U2 from Ui and Uq, % from U2, Ui, and Uq, etc. The first few coeffi- 
cients of the expansion (105) for the eigenvalues u=Ui(a) can thus be obtained from the 
relations : 
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26. Determination of the Eigenvalues Themselves 



We now have to pass from the quantities u=Ui{a), determined by (105), to the parameter 
I which eigenvalue constitutes the complex order of the Legendre function entering in the 
rigorous mode expansion (24). The connection between Ui{a) and I is obtained by combining: 



(a) a relation resulting from (9), for r=a, and (14), viz: 

/(/+1) 



kla' 
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(108) 



(b) the expansion (78) for r=a, which reads: 
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In the latter expansion we substitute for each power of {ut{a)lkoa}^'^ the proper power series 
following from (105) with the aid of a binomial series. We thus arrive at: 
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The quantity Cf is obtained next with the aid of (108), while its square root Ci'^l/(k(,a) 
can be expanded into powers of {koa)~^'^ by a further application of the binomial power series. 
The final result starts as follows: 
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The convergence of this expansion can be proved in special cases, ojie of \\hich will be treated 
iji part III. 

According to (40) the exponentials exp (ild) determine the beliavior of the individual modes 
at large distances. In particular the factor 



, / lod I 



p-Im IqB 



is decisive for the attenuation at these distances, if Iq represejits th(^ (ugejivalue witli the smallest 
imaginary part. Obviously, the eigenvalues I depend on tlie Jiumerical coefficients jj and the 
quantities Uj. In view of (107) the latter follow in succession from tlie other Jiumerical coeffi- 
cients and from Vq. In its turn Uq has to be determined from (104) which (equation, to be 
discussed now iii more detail, represents the zero-order approximation for large values of 

27. Equation for the Zero-Order Approximation in the Case of the Electrical 

Solution 



The equation (104) can be put hi the form: 
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The quotient V JBqq can be reduced as follows to the physical data. Accordhig to section 
6 the parameter V' = Vl{kQay^ is given for the electric solution as follows in terms of the atmos- 
pheric refractive index n{a), its derivative n\a) and the refractive index rie of the earth: 



v'=i{h<iy'Ma)^^^- 



an' (a) 



ni {hay^ia) 
In the second term we can introduce the definition of the effective earth's radius ^eff, viz, 

1 _1 , n'{a) 
Hence, also : 
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47 



1/3 



4 



{hay 



(110a) 



(111) 



Further, we find with the aid of (91), (89), and (88): 

/O \l/3 

<2^') • (112) 
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On the other hand, Ai can be expressed, with the aid of (17), (15), and (6), in terms of higher- 
order derivatives of the profile function n{r) at the earth^s surface. After some reduction 
we find (again for the electric solution): 



Aj — - 



i{(-«S(?--Q},.; 



and a corresponding value for ^oo- 
In terms of the new constant: 

^eff=3T7^' (114) 

which can thus be computed with the aid of (111), (112), and (113), the equation (110) reads: 

TTCl)/ \ .27r/3 
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Substituting Uo={—2Ty^^J3, we get the alternative form: 



eff 



V-2r b. ..,„( =- V' (115) 
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which is identical with a basic equation in the conventional diffraction theory for a homo- 
geneous atmosphere provided that d is replaced there by d^s. [See H. Bremmer, Terrestrial 
Radio Waves, page 43, eq (21a).] The eigenvalues ti satisfying this equation are in the first 
quadrant, and —2ri is defined as a quantity in the —2^^ quadrant. 

The techniques for solving (115) in the case of a homogeneous atmosphere can now also 
be applied here. In particular, the following two expansions could be derived from a Riccati 
equation equivalent to (115) [loc. cit., page 45, formulas (27) and (28)]: 
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The numbers r^.o and r^.^o are defined such that Xi^o=—~ {—2ti^^^'^ and X;,«, =— - (— 27^,^0)^^" 

represent the zeros of Jin{x)-\-J-iii{x) and of J2/3(x)—J-2/3(x). We can obtain the corre- 
sponding expansions for (— 2t)^'^ and thus for u^. The resulting expressions, viz, 
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(116) 
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are particularly useful for computations concerned with small and large values of |5eff|, 
respectively. 

28. Complete Evaluation of the Eigenvalues C^ 



The final complete procedure for arriving at the numerical values of the parameters Ci 
in the case of a vertical dipole can be summarized as follows. First of all, the physical data 
to be given are: 

(a) the function n{r) representing the smooth profile of the refractive index of the atmos- 

phere, 

(b) the refractive index rie of the earth. It is given as follows in terms of the dielectric 

constant e of the earth relative to the atmosphere, and the earth's conductivity cr in 
electrostatic units: 



CO 



The complex parameter V then follows from (111), the effective earth radius being- defined 
by (110a). The parameter Ai can be evaluated next with the aid of (113). It leads to a 
corresponding value of Bqq according to (112), whereupon (114) yields the parameter b^u. 
The successive eigenvalues Uqj can be computed next with the aid of (116). This concludes 
the determination of the zero-order approximation (extended W.K.B. approximation). 

The further corrections depend on a sequence of real parameters startmg with the quan- 
tities A^;. With the aid of (17), (15), and (6) we verify that the latter can be represented by: 



in which 



A.=«^{|.M^w },.^5«"^¥^{^^™5©}.; ^'''^ 



M{r)=-n{r) 
a 



is the modified refractive index of the atmosphere. The second term of (117) is small compared 
to the first one for most relevant frequencies, but it may become important at the lower end 
of the spectrum. 

With the aid of (83) we derive from A^,- the new parameters ft, next with the aid of (81) 
from ft the quantities 7^-. The quantities 7^- and A^ now being known, we derive with the aid of 
(85) the p/s, whereas the g^^'s follow from (88) if we just know the A^'s. Having derived both 
the qnjS and ;^>^-'s, the #„/s can be computed with the aid of (89). The knowledge of the tn/^ 
and g„/s enables the derivation of the Bt^s from (91). Next we compute the Cn/s with the 
aid of (95). The a^/s can then be computed according to (102). AU numerical coefficients 
necessary for the derivation in succession of Ui, U2^ u^ • • • from (107) are then known. 
Finally, Ci can then be computed with the aid of (109), remembering the definition (6) for n^u* 

29. Limiting Cases r'=0 and r' = oo 

These special cases are of particular interest. First of all, the evaluation of the final 
expansion (109) becomes much simpler here than for some finite value of V\ On the other 
hand, the eigenvalues are generally expected to be situated in tlie complex plane between 
those corresponding to the two limiting cases. 

The value r' = is approached for a perfectly conducting earth. In fact, the first term of 
(HI) then vanishes since 71^= 00 while the second term is small but for extreme low frequencies. 
The other limiting case r' = oo is characteristic for short waves for which the earth behaves 
mainly as a dielectric; the factor n(a) (rz./— l)V^e^ in the first term of F' then is of the order of 
unity so as to have a modulus of this term which is of the order of the very large quantity 
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We first deal with the limiting case r'=0. The relations (107) representing ui, U2. 
terms of Uo here reduce to: 
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We substitute these expressions into the eigenvalue expansion (109) wliile approximating 
rieffia) by unity. This approximation amounts to an error which is certainly smaller than that 
introduced by taking r' = for some practical conditions. The resulting expansion reads: 
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We next pass to the other limiting case r' = oo which is simpler. Evaluating (107) we 
first find 1^1 =0. However, when determining the limit for F'— >oo in the successive terms for 
the higher order u/s, we have to take into account the approximation 
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for large V\ We then obtain the following expressions : 
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Substituting these values into (109), again taking neff((i) = l, the eigenvalue expansion for 
r' = oo proves to start with: 
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Examples of the limiting eigenvalue expansions (118) and (119) will be worked out nu- 
merically in part III. The first few terms of the expansions (118) and (119) are rather simple 
even when all its parameters are expressed explicitly in terms of the profile constants A^. The 
c()in|)lete evaluation of these parameters results in expansions starting as follows: 
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30. Plane Stratification as a Limiting Case 



The theory for a concentrically symmetric medium can be reduced to that for a medium 
with a flat stratification (refractive index dependent on an orthogonal coordinate z instead of 
the radial distance r). The determination of solutions and eigenvalues of the ''flat height-gain 
differential equation' ' has been investigated by various authors who discuss extensions of the 
W.K.B. approximation, comparable to those of the present paper. We mention in particular 
articles by Cherry [1950], Imai [1956], and Erdelyi [I960]. Recently the problem has been 
approached in detail by Logan [1961] for the eigenvalues. However, we only considered the 
spherical problem in view of its direct applicability to far-distance troposplieric propagation. 
As a matter of fact, results for a flat stratiflcation may be derived from those given liere with the 
aid of a limiting procedure. This will be explained now. 

Our basic differential equation has the following form, in view of (8), (9), (15), and the 
coefficients (17) of the Taylor expansion for ilffff(r) around r=a: 



We substitute 



r=a+z; Aj=a%; rJ,{r) = F{z); 1(1 + 1) ^P=\'a', 



(122) 
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and next take the limit for a— >oo. The differential equation resulting reads: 

On the other hand, the boundary condition (23) becomes, with T=ya, 

g'=yFat2=0. (124) 

In view of the relation T=(kiiay^ V, and the parameter 

Boo={^Uo) =(^2^'V =V2N =12^^-; ' 
the equation (104) defining the roots Uq may now be represented b}^-: 
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The solution of the differential equation (123), with the boundary condition (124) at 2=0 
and the radiation condition at 2==+oo^ constitutes a very general eigenvalue problem. Its 
properties can be derived from those in the present paper by substituting (122) into the latter, 
and then passing to the limit for a-> oo , As an example we can give the expansions which 
correspond in the spherical case to (120) and (121): 
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